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ABSTRACT 

cn 

^— ' Although tidally-locked habitable planets orbiting nearby M-dwarf stars are 

pq among the best astronomical targets to search for extrasolar life, they may also 

?H be deficient in volatiles and water. Climate models for this class of planets show 

^ atmospheric transport of water from the dayside to the nightside, where it is 

■ precipitated as snow and trapped as ice. Since ice only slowly flows back to 

(NJ the dayside upon accumulation, the resulting hydrological cycle can trap a large 

amount of water in the form of nightside ice. Using ice sheet dynamical and ther- 

modynamical constraints, I illustrate how planets with less than about a quarter 

the Earth's oceans could trap most of their surface water on the nightside. This 

p I would leave their dayside, where habitable conditions are met, potentially dry. 

Q The amount and distribution of residual liquid water on the dayside depend on a 

-!—> variety of geophysical factors, including the efficiency of rock weathering at reg- 

c3 ulating atmospheric CO2 as dayside ocean basins dry- up. Water-trapped worlds 

with dry daysides may offer similar advantages as land planets for habitability, by 

►^ contrast with worlds where more abundant water freely flows around the globe. 

(N 

.• 1. Introduction 

■^ 

O 

CO Tidally-locked terrestrial planets orbiting in the habitable zone of nearby M-dwarf stars 

are among the best astronomical targets for discovery and atmospheric characterization (e.g., 

Charbonneau & Deming 2007; Montgomery & Laughlin 2009; Seager & Deming 2010; Rojas- 

Ayala et al. 2013). Recent estimates based on Kepler results are encouraging in suggesting 

a high frequency of such planets around M-dwarfs (Dressing & Charbonneau 2013; Fressin 

et al. 2013; Kopparapu 2013; Morton & Swift 2013). 

The climate on this class of planets is peculiar and interesting because of the perma- 
nent day-night insolation pattern many are expected to experience by the end of their tidal 
evolution, once they are captured in 1:1 spin-orbit resonance (e.g. Rodriguez et al. 2012; 
see Correia & Laskar 2011 for a review). Provided that large enough surface pressures are 
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present, nightside atmospheric collapse can be avoided (Joshi et al. 1997). Various other 
studies have illustrated how habitable planets around M-dwarfs will typically find themselves 
in a climate regime with hot daysides, cold frozen nightsides and an atmospheric circulation 
connecting the two hemispheres (Joshi 2003; Merlis & Schneider 2010; Heng et al. 2011; Kite 
et al. 2011; Pierrehumbert 2011; Wordsworth et al. 2011; Edson et al. 2011,2012; Leconte 
et al. 2013). 

Two aspects of the climate problem on these planets which have received less attention, 
however, are the possibilities that they are deficient in volatiles and that a significant fraction 
of their surface water inventory could be trapped as ice on their nightside. Indeed, Lissauer 
(2007) and Raymond et al. (2007) have argued on the basis of standard terrestrial planet 
formation scenarios that volatile deficiency may be a key attribute of these planets (see also 
Ogira & Ida 2009 for an alternate view). Even planets formed with abundant water could 
lose much of their surface inventory via atmospheric escape during an intense tidal heating 
phase preceding their spin-orbit resonant capture (Barnes et al. 2013). Heath et al. (1999) 
and Joshi (2003) have considered the nightside ice trap issue but they mostly emphasized 
how the ice layer would experience melting at its base on worlds with abundant water, thus 
making the ice trap inefficient. Recently, Leconte et al. (2013) have also discussed the 
nightside ice trap but mostly in the specific context of the climate bi-stability of planets 
found near the inner edge of their habitable zone. 

Here, I interpret volatile deficiency as likely implying a small surface water inventory 
for these planets and I argue that the ice trap on the nightside of water-deficient planets 
around M-dwarfs could in fact be efficient and have a profound effect on their climate. In 
^ I present new simulations of the climate on habitable planets around M-dwarfs which 
illustrate the rapid development of a nightside ice trap for surface water. In ^ I develop the 
argument about the effective nature of this ice trap on sufficiently water-deficient planets, 
before discussing some implications for their climate and concluding in ^ 



2. Climate Simulations and Ice Trapping 

I use PlanetSimulator, a fiexible Earth-System simulator of intermediate complexity 
developed at the University of Hamburg]^ to study the climate of tidally-locked habitable 
planets around M-dwarfs. PlanetSimulator relies on a pseudo-spectral atmospheric dynami- 
cal solver coupled to an accurate radiative transfer scheme. It includes detailed formulations 
to simulate the hydrological cycle on Earth-like planets, such as water evaporation and pre- 
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cipitation, diagnostic cloud formation, a SOm slab model for the ocean and a thermodynamic 
sea ice model (see Fraedrich et al. 2005 for details). Although these model elements have 
been tuned to simulate the Earth's climate, they are expected to remain reasonably accu- 
rate for climates deviating only modestly from Earth's conditions, as is the case here. I use 
PlanetSimulator in its T42L10 formulation (128 longitude nodes, 64 latitude nodes and 10 
vertical levels). 

I modified PlanetSimulator to enforce a permanent day-night insolation pattern and I set 
the planetary obliquity and orbital eccentricity to zero in all models, for consistency with the 
tidally-locked hypothesis investigated here. I adopt an Earth-like insolation of 1365 W m~^ 
in all models and a planet rotation period Prot = 240 h by default, which corresponds to 
a tidally-locked planet around a ~ 0.3 Mq M-dwarf (Edson et al. 2012). I also removed 
orography (relief over continents) from all models to make the results more general and after 
tests showed that it does not have a major impact on the main results. I use the current 
distribution of Earth continents in models with continents. All other model parameters take 
their standard Earth values, unless otherwise statedjj 

Table [l] shows the list of models I have explored in this work. All models are run for 40 
Earth months, which is sufficient for the surface temperature field to reach a steady-state 
under permanent insolation, starting from Earth-like initial conditions. Earthl, Earth2 and 
Earth3 are models with standard Earth parameters, the default rotation period Pj-ot = 240 h 
and a permanent insolation pattern centered on longitude deg (Greenwhich), +90 deg 
(Bangladesh) and +180 deg (central pacific ocean), respectively. Earth. lowP is similar to 
Earthl except that the total atmospheric surface pressure is reduced to 0.3 barj^ With 
a constant CO2 fraction of 360 ppm, this also corresponds to a CO2 atmospheric content 
reduced by about a third. Earth. lowC02 and Earth. highC02 are similar to Earthl except 
for CO2 fraction of 36 ppm and 3600 ppm, respectively]^ Earth. fast and Earth. slow are 
similar to Earthl except for their rotation periods Prot = 120 h and 480 h, respectively. 



^In particular, I have not accounted for the reduced ice and snow albedos expected on planets orbiting 
M-dwarfs (Joshi & Haberle 2012). This is expected to have only a minor quantitative impact on the main 
results of interest here because only a small fraction of the icy surface is actually exposed to insolation (see 
also Leconte et al. 2013). 

^Volatile deficiency could in principle also result in a reduce N2 atmospheric content, which is captured 
by the lower surface pressure in this model. Note that improved modeling of low surface pressure scenarios 
would entail verifying the accuracy of the PlanetSimulator radiative scheme under these conditions, which 
was not done here. 

""^Lucarini et al. (2013) suggest that PlanetSimulator's radiative transfer scheme should remain reasonably 
accurate over this range of CO2 atmospheric content values. 



SuperEarth is similar to Earthl except that it uses larger values of the surface gravity and 
the planetary radius, corresponding to a super-Earth with a mass ~ 10 M^ (see Table [I] 
Valencia et al. 2006). Finally, Aquaplanet is similar to the Earthl-3 models except that all 
the continents have been replaced by a global ocean. 

I find that the general circulation regime in all these models is qualitatively consistent 
with what has been reported in the literature before for other tidally-locked habitable planet 
models (e.g., Joshi et al. 1997; Merlis & Schneider 2010; Heng et al. 2011; Edson et al. 
2012). The majority of models listed in Table [I] exhibit a zonally-averaged zonal wind 
profile comparable to that of Earth, with two eastward jets streams at altitude and negative 
wind speeds at the equatorial ground. In three models, however (Earth. lowP, Earth. slow 
and SuperEarth), the general circulation is in a sup errot at inaj regime with positive zonal 
winds at the equatorial ground and a rather broad eastward jet at altitude. Maximum 
zonally-averaged wind speeds range from ~ 17 m s^^ (SuperEarth model) to ~ 60 m s^^ 
(Earth. highC02). Interestingly, these results are not entirely consistent with the study of 
Edson et al. (2012), who reported faster wind speeds in their superrotating models and 
dynamical transitions to the superrotating regime at Prot — 100-101 h for land planets and 
Prot — 72-96 h for aquaplanets (i.e., at faster rotation rates than found here). It would be 
interesting to understand the origin of this discrepancy better. I simply note for now that 
it does not appear to have a major impact on the main results of interest here for nightside 
ice since all models in Table [T] exhibit comparable properties in that respect. 

Figure [l] shows representative end results from the Aquaplanet model. The surface 
temperature field shows that liquid water conditions are confined to the dayside of the planet 
(see also Figure [i]), making this model climate comparable to the eyeball climate discussed 
by Pierrehumbert (2011). On the nightside, the coldest surface conditions are offset from 
the poles and the equator is the hottest, as the region where the near-surface wind is the 
most efficient at carrying heat from the dayside and coupling it with the nightside surface. 

Precipitable water in the atmosphere is also largely confined to the dayside, where evap- 
oration strongly peaks at the substellar point. The circulation is able to advect, preferentially 
eastward, a modest amount of moisture to the nightside, where it is precipitated as snow. 
Most of the snow precipitation occurs in an annulus on the dayside, but the weak snow pre- 
cipitation rates on the nightside are nevertheless interesting in terms of the ice accumulation 
rates they imply (see below). As noted before by Merlis & Schneider (2010) in comparable 



^Superrotation refers to a regime with positive equatorial winds, which imply "eddy" momentum transport 
towards the equator beyond what any type of axisymmetric meridional circulation can achieve (e.g., Showman 
et al. 2011). 



simulations, strong net precipitation occurs in the vicinity of the substellar point, with net 
evaporation over most of the dayside. All of the properties illustrated in Figure [T] are qual- 
itatively shared by the other models listed on Table [T| which suggest that they are rather 
general for tidally-locked habitable planets around M-dwarf stars. 

For each of the models. Table [T] also lists the snow precipitation rate averaged over the 
planetary nightside. Values range from 0.3 to 0.01 mm day~^, in broad agreement with 
the 0.11 mm day~^ nightside rate reported by Joshi (2003) Jj This range largely reflects 
the efficiency of the circulation at advecting moisture to the nightside, which is noticeably 
weaker in the Earth. slow model (see also Merlis & Schneider 2010). Nevertheless, with 
such precipitation rates, a layer of 1 km of ice would take between 3 x 10^ and 10® yr, and 
typically ~ 10^ yr, to accumulate on the nightside (accounting for the factor ~ 3 difference 
between the density of snow and ice). These timescales are quite short and they suggest 
a rapid change in the hydrological cycle of planets with no more than a few km of ocean 
(such as Earth), following their capture into 1:1 spin-orbit resonance. It is thus important to 
understand the nature of this ice layer in detail if one is to address the climate of habitable 
planets around M-dwarfs. 



3. Nightside Ice Sheet 

3.1. General Considerations 

Figure [2] illustrates schematically three possible configurations for the nightside ice layer 
on a tidally-locked planet. For the sake of argument, a uniformly fiat ocean basin has 
been assumed and continents have been ignored, though the consequences of lifting these 
assumptions are reconsidered in ^ 



As discussed quantitatively in §3.3[ a conducting ice layer carrying a given geothermal 
heat fiux has a limited thickness because melting eventually occurs at its basej^ This limiting 
thickness is smaller for warmer surface temperatures at the top of the ice layer. Tidally- 
locked planets with relatively mild nightside temperatures and/or a large surface water 
inventory will thus have their nightside ice layer fioat over a sub-glacial ocean connecting to 
the dayside ocean. This scenario, which is the one emphasized by Health et al. (1999; see 



^Note that Joshi's simulations are not fully consistent for M-dwarf habitable planets in that they assume 
an Earth-like planetary rotation rate (Piot = 24 h). 

^At least for ice pressures ^ 100s MPa, which is the relevant regime for the water-deficient planets of 
interest here. 



also Joshi 2003), is depicted in the first "ice shelf configuration" shown in Fig. M In this 
configuration, rain and snow precipitate on the nightside, as a result of moisture advection 
from the dayside (gray arrows). Snowfall does not add to the ice layer thickness, however, 
because the thickness is limited by basal melting. The steady-state hydrological cycle is 
closed with a return fiow of water to the dayside via the sub-glacial ocean (yellow arrows). 
The ice-|-ocean layer thickness must be essentially uniform over the entire globe to satisfy 
lateral pressure equilibrium at the base, or rapid water fiows would ensue to enforce this 
equilibrium. 

As depicted by the second "ice sheet configuration" in Fig.|2} on a planet with less surface 
water and/or colder nightside temperatures, it becomes possible for part of the ice layer not 
to exceed the critical thickness for melting. This region of the ice layer is then grounded, as 
shown in between the two radial dotted lines in Fig. [2j Rain and snow precipitation act the 
same as in the previous configuration, with the important exception of the grounded ice sheet. 
Ice is able to accumulate in this ice sheet and sustain a sizable pressure gradient at its base 
(indicated by a bulged ice layer in Fig. |2]) because of the considerable ice viscosity. Ice fiows 
down the basal pressure gradient in the ice sheet and, in steady state, the hydrological cycle 
must satisfy a new mass balance: over the grounded ice sheet, ice accumulation from snowfall 
must balance the fiow of ice exiting the ice sheet at its boundary with the sub-glacial ocean. 
Since ice accumulation is made possible by the enormous difference of viscosity between the 
grounded ice flow (in the ice sheet) and the water flow (in the ocean) covering the rest of 
the globe, an ice sheet model describing the viscous flow of ice is necessary to evaluate the 
amount of ice present in the ice sheet. 

The third "water-trapped conflguration," which is of particular interest to the present 
study, is essentially an extreme version of the "ice shelf conflguration." On a planet with even 
less surface water and/or even colder nightside temperatures, one can in principle have the 
great majority of the surface water trapped in the slowly flowing ice sheet, with little water 
present in liquid form on the dayside. A flnite liquid dayside inventory remains a necessity 
to close the hydrological cycle. Indeed, it is the evaporation, advection and subsequent 
precipitation of this dayside water on the nightside that must balance the flux of ice flowing 
back to the dayside, but the size of this dayside hquid reservoir could in principle be very 
small, as depicted in Fig. [2j To evaluate how much water can potentially be trapped on the 
nightside in this conflguration, I now turn to quantitative ice sheet models. 



3.2. Isothermal Ice Sheet Models 

Thermo-mechanical models of ice sheets exist (e.g. Bueler et al. 2007; Fowler 2011) 
but I find it convenient for the present purpose to consider separately mechanical and ther- 
modynamical constraints to gain insight into the behavior of nightide ice on a tidally- locked 
planet. The use of a simple, isothermal ice sheet model can be justified by noting that a 
comparable diffusion equation can be derived for the non-isothermal case (e.g., Fowler 2011) 
and that ice sheets generally exhibit strong vertical shear preferentially near their bottom 
boundary, i.e. over a limited range of temperatures. The evolutionary equation satisfied by 
a thin isothermal ice sheet with a fiat base and a no slip boundary condition at its bottom 
is (Fowler 2011) 
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where h is the ice thickness, p = 917 kg m~'^ is the mean density of ice, g is the surface 
gravity, Aq = 10~^® Pa yr~^ is a typical value of the fiow rate constant (related to the ice 
viscosity near the bottom of the ice sheet), n = 3 is the typical value used in Glen's ice 
fiow law to model the non-Newtonian nature of the ice fiuid and a is the ice accumulation 
rate (assumed constant to model uniform snowfall over the ice sheet). This equation is a 
non-linear diffusion equation for the ice thickness, h, satisfying the conservation of mass and 
momentum in the shallow ice approximation {h much less than the horizontal length-scale of 
the ice sheet). The quantity in bracket is the ice flux, with an effective diffusion coefficient 
that depends strongly on the ice thickness and its horizontal gradient. 

Equation ([l]) admits analytical solutions (e.g., Bueler 2003) and I shall use such a simple 
solution here. I solve the adimensional, steady-state axisymmetric version of Eq ([I]) for the 
standard n = 3 value in Glen's law. 
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where r = R/L is the adimensional cylindrical radius from the ice shelf center (located 
in i? = 0), r = 1 marks the outer edge of the axisymmetric ice sheet (defined by i? = L 
dimensionally) , s = h/Z is the adimensional ice thickness, and 
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is the ice thickness scale factor. 

I integrate this equation a first time with respect to r and use the central boundary 
condition ds/dr = in r = (the ice "ridge"). I integrate a second time with respect to r 
and use a finite, arbitrary thickness /lout for the boundary condition in r = 1. This boundary 
condition guarantees a finite ice fiux at the point where the ice ceases to be grounded and 
reaches fiotation over water (indicated by the radial dotted lines in the middle panel of 
Fig. pi), at which point the ice sheet equation ceases to be valid (fowler 2011). The resulting 
dimensional solution for the ice sheet thickness is 
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This solution generalizes the singular case /lout = already known in the literature (e.g., 
Bueller 2003). The case /lout > is more satisfactory for the present application in that it 
permits a finite ice fiux at the outer edge (where fiotation is reached) to balance the ice 
accumulation rate integrated over the entire ice sheet. 

Figure IS] shows a few ice thickness profiles derived from Eq. (|4]). The solid line is 
the solution for an ice sheet of L = 1000 km in radius, a uniform ice accumulation ratqj 
a = 0.1 mm day^^ and a boundary thickness /lout = 1 km. The dash-dotted line corresponds 
to a similar solution except for /lout = 500 m, while the dotted line corresponds to a reduced 
accumulation rate a = 0.01 mm day^^. The dashed line is the solution for an ice sheet of 
L = 500 km in radius (with accumulation rate a = 0.1 mm day~^ and /lout = 500 m). The 
corresponding values of the ice thickness scale factor in these four solutions are Z = 2446, 
2446, 1835 and 1730 m, respectively. 

An exploration of the parameter space of these simple ice sheet solutions shows that 
the value of hout has only a small impact on the global ice profile, as long as /lout ^ Z. The 
scale factor Z itself is a good indicator of the overall ice thickness across the profile (and 
thus the total mass ~ ZL^ in the ice sheet), as illustrated in Fig. ^ Stronger (weaker) ice 
accumulation rates lead to thicker /steeper (thinner/shallower) ice profiles, as expected from 
the ice sheet adopting a steady-state thickness profile such that the local ice fiux at any 
radius carries an amount of ice equal to the integrated accumulation rate within that radius 
(Eq. (|2])). It is this mass balance that sets the typical ice thickness in these models. 



®Note that ice accumulates at a rate approximately 3 times slower than snow precipitates, in proportion 
to the ratio of ice to snow density. 



For ice accumulation rates representatives of the average nightside precipitation in the 
chmate models discussed in ^{ 0.01-0.3 mm/day of snow, see Table [I]), the typical ice thick- 
ness is thus Z ~ 1600-2450 m in an ice sheet of radius L = 1000 km (Eq. p])). The climate 
models suggest the possibility of an ice sheet extending globally over the entire nightside, 
with a radial extent L ~ 10000 km on an Earth-size planetjj For such a large horizontal 
extent, the typical ice thickness reaches Z ~ 4950-7750 m, according to Eq. ([s]). These 
depths would correspond to a considerable amount of ice being trapped on the nightside. As 
we shall now see, however, such a massive ice trap is not expected on the basis of additional 
thermodynamical constraints. 



3.3. Thermodynamic Melting Constraint 

Ice thickness is also limited by the melting condition for ice at high pressure, at the 
base of the ice sheet. In its simplest form, this process is modeled as a vertical, steady-state 
conduction problem for the ice layer transporting a specified geothermal flux to the surface. 
Here, I adopt for concreteness the formalism of Abbot & Switzer (2011), which accounts 
for the important variation with temperature of the ice thermal conductivity. The resulting 
exponential temperature profile for the steady-state conductive ice layer yields a maximum 
thickness before melting, 

'icond = T^ In I 7^ ) , (5) 



-^ V -^ surf , 

where the water ice conductivity is given by k{T) = AT"^, A = 651 W m~^, T is 
the ice temperature in K, F is the thermal flux being conducted (~ 0.09 W m~^ for Earth's 
geothermal flux) j^Tmoit — 260 K is a representative high-pressure water melting temperature 
and Tsurf is the temperature boundary condition at the top of the ice layer. Tsnri is determined 
by surface-atmosphere exchanges which are modeled explicitly in the climate simulations 
described in [|2} 

It is therefore straightforward to translate a surface temperature map into a map of 
maximum ice thickness via Eq. ([5]). An example is shown in Fig. Ill for the Aquaplanet 



^AU else being equal, a super-Earth has ice sheets g ^/^ times thinner from stronger gravity, although a 
global ice sheet could also be thicker from the extra L^/^ horizontal scale factor in Eq. tin. 

^"The geothermal flux may scale up with planet mass to the power one half, so that ice layers would be 
comparatively thinner on super-Earths (Abbot & Switzer 2011). 
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model, assuming Earth's geothermal flux for F . This map shows that the melting constraint 
on the ice thickness is stringent. Most of the nightside ice has a maximum depth ~ 1-1.5 km 
before melting occurs, with a colder region of a few thousand km in extent allowing for thicker 
ice (up to 4.3 km depth). Since most of the nightside exhibits surface temperatures which do 
not permit ice thickness in excess of 1.5-2 km thermodynamically, the considerably larger 
values of ice thickness (5-8 km) that would otherwise accumulate according to steady-state 
ice sheet models cannot be realized in practice: basal melting will prevent such large ice 
buildup. One concludes that, on global planetary scales and for the precipitation rates listed 
in Table [l} the thermodynamic constraint will be the one limiting the amount of ice that can 
be trapped on the nightside of a tidally-locked planet. The actual ice thickness may be close 
to marginal basal melting since steady-state mass balance, as described by simple isothermal 
ice sheet models, would systematically have more ice accumulate, until basal melting limits 
this accumulation p] In what follows, I will adopt this point of view, by using the maximum 
ice thickness described by Eq. ([s]) as a reasonable proxy for the steady-state thickness of ice 
that is actually trapped on the nightside. 

After integration, the ice distribution shown in Fig. |4] amounts to an average depth of 



560 m spread over the entire globe. Similar averages^ for the other climates models listed in 
Table [T] range from 320 to 770 m. This range of values largely reflects how cold the nightside 
is in each model, with the lowest efficiency of heat advection to the nightside found in the 
model Earth. low (with a low atmospheric mass). By comparison, the amount of surface 
water on Earth (essentially in the oceans) corresponds to an average depth of ~ 2700 m once 
spread over the globe (e.g., Charette & Smith 2010). One concludes that the climate models 
listed in Table [T] can trap from a seventh up to almost a third of the Earth's surface water 
on their nightsides, with the majority of models trapping about a quarter to a fifth. This is 
considerably more that the amount of water in ice form on Earth and it is a large enough 
amount that water-deficient planets may be at risk of trapping most of their surface water 
on their nightside. 

It is difficult, however, to make robust predictions for the total amount of water present 
on terrestrial planets on the basis of planet formation models. Simulations of in situ forma- 
tion by Raymond et al. (2007) suggest that habitable zone planets around M-dwarfs could 
form with as little as two orders of magnitude less bulk water than Earth, which would likely 
result in significantly reduced amounts of surface water as well. On the other hand, Ogira 



^^An ice sheet with basal melting will flow more rapidly through sliding. A thermo-mechanical model 
would be needed to capture faithfully this behavior (Fowler 2011). 

^^For simplicity, I do not differentiate between land and ocean when making such averages. The role of 
continents is reconsidered in f|4l 
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& Ida (2009) argue that water could be abundantly present on terrestrial planets in the 
habitable zone of M-dwarfs if planetary migration efficiently brings in bodies from beyond 
the ice line. Depending on details of the protoplanetary disk structure and how it influences 
the planetary migration process, both scenarios could be at work in nature. M-dwarf planets 
could also lose much of their initial surface water inventory via atmospheric escape during 
the intense tidal heating phase that may precede their spin-orbit resonant capture (Barnes 
et al. 2013). I will adopt the simple view here that at least some terrestrial planets in the 
habitable zone of M-dwarfs could be quite deficient in surface water. 

Therefore, among the population of habitable Earth-like planets around M-dwarfs, those 
who possess only a fraction of the Earth's surface water inventory could trap the great 
majority of their water in the form of nightside ice. Although the dayside of such water- 
trapped worlds may end up being quite dry, it is worth emphasizing again that a closed 
hydrological cycle requires liquid water to remain present on the dayside to compensate for 
the return flow of ice from the nightside (see Fig |2] and discussion below). Evaluating exactly 
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how much liquid water is present on the dayside of a water-trapped world is challenging 
and, as I shall briefly describe below, it is closely related to other long-term climate issues 
for such planets. 



4. Discussion and Conclusion 

Using a combination of climate and ice sheet models, I made the case for tidally-locked 
habitable planets around M-dwarfs being able to trap most of their surface water in the 
form of nightside ice, provided they are sufficiently volatile-deflcient. Quantitatively, planets 
subjected to an Earth-equivalent insolation and possessing less than about a quarter the 
Earth's surface water inventory could rapidly flnd themselves into a water-trapped conflgu- 
ration, following their capture into 1:1 spin-orbit resonance. The critical amount of surface 
water below which a water-trapped conflguration can occur will probably vary with the 
magnitude of insolation received by such planets, although this issue cannot be easily disen- 
tangled from the atmospheric CO2 content and the climate being regulated by global-scale 
weathering, which has a potentially strong dependence on the continental distribution for 
tidally-locked planets (Edson et al. 2011). 

Clearly, the main argument constructed in this work relies on a number of assumptions 
and approximations made in the climate and ice sheet models, which all deserve further 



^■^For example, Leconte et al. (2013) have presented an idealized model for a water reservoir present at 
the edge of the icy region, without accounting for constraints on the nightside ice flow. 
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scrutiny. For example, oceanic heat transport to the nightside, which is not accounted 
for in either type of models, could provide an extra heat source that may act to limit the 
ice thickness relative to the results discussed above, ^ In the limit where this transport 
is efficient, significant ice accumulation may only occur on the fraction of the planetary 



nightside that is covered by continents^ thus reducing the size of the ice trap. This is an 
example of how the detailed topography and ocean basin configuration on a water-trapped 
world can have a strong impact on its climate. 

Another example of the important role played by the continental distribution and ocean 
basin configuration concerns the amount of residual liquid water left on the dayside of a 
water-trapped world. In the climate models discussed in ^ the ability to precipitate snow 
at the rates listed in Table [T] largely relies on evaporation over oceans and moisture transport 
to the nightside by the atmosphere. In a water-trapped configuration, as the dayside oceans 
gradually dry up, their evaporative surface will shrink in a way that depends on the detailed 
ocean basin configuration. The atmospheric moisture transport to the nightside may also 
be affected when a significant lowering of the dayside sea level occurs. If this were to result 
in substantially decreased nightside precipitation rates, the nightside ice sheet would tend 
to accumulate less ice than evaluated in ^ (see, e.g., scaling with a in Eq. (|3])). Eventually, 
the ice layer depth would not be limited by melting at the base but simply by mass balance 
at very small accumulation rates. This regime, which would require snow precipitation rates 
orders of magnitudes smaller than the values listed in Table [T| may happen on some water- 
trapped worlds but not others, when dayside oceans dry up. Regardless, this reasoning 
illustrates how a reliable evaluation of the amount of liquid water left on the dayside of a 
water-trapped world must rely on a closed hydrological cycle, i.e. one that accounts for two- 
way exchanges between the dayside and the nightside water reservoirs, in a manner more 
explicit than the isolated ice sheet models considered in ^ 

The amount of residual liquid water on the dayside of water-trapped worlds will also 
depend on the efficiency of rock weathering at regulating the CO2 content of the atmosphere. 
Abe et al. (2011) argue that weathering will continue to operate on land planets but the 
efficiency of the process may differ from what it is on Earth, depending on the availability of 



^^Under particular conditions, tidal dissipation could contribute an additional heat flux acting to limit the 
nightside ice layer thickness (e.g., Barnes et al. 2013; Leconte et al. 2013). 

^^I note that the literature on Earth's True Polar Wander suggests that continental distribution and 
planetary moments of inertia are not strongly coupled attributes (e.g., Kirschvink et al. 1997, Raub et. 
al. 2007). This also suggests a lack of correlation between the continental distribution and the day/night 
hemispheres that end up being locked into 1:1 spin-orbit resonance by tides on a close-in planet. 
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surface water for different land regions, ^^ On a planet with dried-up dayside ocean basins, 
more rock will be exposed to the atmosphere but the opportunity for ocean weathering would 
also largely disappear, with unclear consequences for the global weathering rate (e.g.. Abbot, 
Cowan & Ciesla 2012). Quantitatively, the models Earth. lowC02 and Earth. highC02 listed 
in Table ?? exhibit ~ 150 m differences in ice their equivalent thicknesses, relative to the 
model Earthl with 360ppm of CO2. This makes it clear that even modest adjustments to 
the atmospheric CO2 content by weathering can have a significant impact on the residual 
dayside surface water inventory of a water-trapped world. 

A related issue is the fate of the large CO2 reservoir stored in the ocean (about 50 
times that in the atmosphere on Earth). Although this reservoir may shrink in proportion 
to the total ocean mass on volatile deficient planets, if one assumes a fixed CO2 ocean 
concentration, it could still amount to a significant and sudden CO2 increased if released 
in the atmosphere when dayside oceans dry-up on a water-trapped world. This excess CO2 
would presumably be weathered on long-enough timescales, so that the long-term climate 
would remain determined by detailed balance between the hydrological cycle and the carbon- 
silicate cycle. 

The trapping of the majority of the surface water inventory in the form of nightside ice 
does not necessarily constitute a challenge to the habitability of water-trapped worlds. Abe 
et al. (2011) found more extended habitability limits for land planets than for aquaplanets 
and they illustrated how land planets are better at holding on to their water, which is lost 
through atmospheric escape, even if they start with less water than aquaplanets (see also 
Leconte et al. 2013; Zsom, Seager & de Wit 2013). Water-trapped worlds with dry daysides 
and frozen nightsides may thus offer similar advantages as land planets for habitability, in 
comparison with worlds where more abundant water is free to flow around the globe (first 
configuration shown in Fig. |2|. A water-trapped configuration may, however, lead to more 
strongly localized habitable regions, depending on the availability and detailed distribution of 
surface water on the dayside. It would thus be interesting to determine to what extent water- 
trapped and water-free configurations can be discriminated from one another observationally. 
This may be challenging if the surface conditions are not directly accessible (e.g., Benneke 
& Seager 2012) and if the stratospheric water content probed by transmission spectroscopy 
is not too different between these two configurations (e.g.. Fig. 12 of Joshi 2003). 

The author is grateful to D. Valencia for many insights on the geophysics of this problem 



^^For example, desertification could occur in regions where evaporation dominates over precipitation, such 
as the dayside annular region away from the substellar point with negative net precipitation in Fig. [T| 
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Model 


Key Parameters 


Nightside Snow 


Equivalent Ice 






Precipitation (mm/day) 


Thickness (m) 


Earthl 


Pj,oj = 240 h, Stellar longitude = deg 


0.12 


470 


Earth2 


Prot = 240 h, Stellar longitude = +90 deg 


0.11 


440 


Earths 


Pj.„t = 240 h, SteUar longitude = +180 deg 


0.26 


400 


Earth.lowP 


Earthl except -Psurf = 0.3 bar 


0.27 


770 


Earth.lowC02 


Earthl except C02 = 36 ppm 


0.10 


600 


Earth.highC02 


Earthl except C02 = 3600 ppm 


0.16 


320 


Earth. fast 


Earthl except P^ot = 120 h 


0.29 


340 


Earth. slow 


Earthl except Prot = 480 h 


0.009 


490 


SuperEarth 


Earthl except g = 29.63 m s"^^ ^^ = hqOO km 


0.07 


400 


Aquaplanet 


Earthl without continents 


0.06 


560 



Table 1: Key parameters of climate simulations and results. 
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Fig. 1. — End snapshots from the Aquaplanet modeL Clockwise, from upper left: sur- 
face temperature (K), precipitable water (kg m~^), net precipitation (precipitation minus 
evaporation in m s~^, with an arbitrarily skewed scale to highlight low values) and snow 
precipitation (m s~^, logarithmic scale with zero reset at 10~^*^). The substellar point is 
located at zero longitude and latitude in each plot. Surface water conditions are confined 
to the dayside and the coldest regions on the nightside are offset from the poles. A modest 
amount of moisture is advected to the nightside, resulting in weak snow precipitation rates. 
Net evaporation occurs on the dayside in an annulus from ~ 20 deg to ~ 50 deg from the 
substellar point, with net precipitation over the remainder of the globe. 
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Fig. 2. — Schematic plot of three ideahzed configurations envisioned for the surface water 
inventory on a tidally- locked terrestrial planet. A view from the pole is shown, with the 
nightside up and the dayside down. Exaggerated thicknesses are adopted for the ocean 
(blue) and ice (white) layers. In the ice shelf configuration, the ice layer builds up on the 
nightside but melting occurs at its base well before the ground is met. In the ice sheet 
configuration, the ice layer is partly grounded over the coldest regions on the nightside 
(between the two radial dotted lines) and it connects to an ice shelf further away. In the 
water-trapped configuration, most of the surface water is trapped in the nightside ice sheet, 
with comparatively little water on the dayside. In all cases, the hydrological cycle is closed 
by advection of atmospheric moisture from the dayside to the nightside, where it is deposited 
as snow (gray arrows), and a return fiow of water to the dayside in the ocean (yellow arrows). 
Planets with less surface water and/or colder nightside temperatures, are more likely to be 
in the water-trapped configuration. 
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Fig. 3. — Axisymmetric solutions for a thin isothermal flowing ice sheet. The solid line is 
the solution for an ice sheet of 1000 km in radius, a uniform ice accumulation rate of 0.1 mm 
day^^ and an outer boundary thickness value of 1 km. The dash-dotted line corresponds to a 
similar solution except for an outer boundary thickness value of 500 m, while the dotted line 
corresponds to a reduced ice accumulation rate of 0.01 mm day^^. The dashed line is the 
solution for an ice sheet of 500 km in radius (with accumulation rate of 0.1 mm day"^ and 
outer boundary thickness value of 500 m). These steady-state ice profiles carry an outward 
fiux of ice at each radius that is equal to the integrated ice accumulation rate within that 
radius. 
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Fig. 4. — A representative map of the ice thickness distribution (in meters) permitted by 
the thermodynamic melting constraint in the Aquaplanet model, for a geothermal flux equal 
to that of Earth. Most of the nightside ice has a maximal depth ~ 1-1.5 km, with a colder 
region of a few thousand km in extent permitting a thicker ice layer (up to 4.3 km). The 
equivalent depth of this ice distribution, averaged over the entire globe, is 560 m. 



